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1. Introduction 



The possibility of superluminal photon propagation in gravitational fields is an intrigu- 
ing prediction of quantum field theory in curved spacetime. It raises fundamental questions 
about the realisation of causality and is potentially of great importance for early-universe 
cosmology. 

In the context considered here, superluminal propagation was originally discussed in 
a pioneering paper by Drummond and Hathrell |1J , who investigated the effects of vacuum 
polarisation on photon propagation in a variety of classical curved spacetimes, including 
Schwarzschild, Robertson- Walker and weak-field gravitational waves. They showed that 
in general it is possible to find directions and polarisations for which the photon velocity 
exceeds the fundamental constant c. This work has subsequently been generalised to other 
examples of background spacetimes, including Reissner-Nordstrom and Kerr black holes 
In ref.Q, we presented some further theoretical analysis, including the formulation 
of a polarisation sum rule and a horizon theorem, based on the original work of [0 . In this 
paper, we extend this development and present results for a phenomenologically impor- 
tant spacetime, the Bondi-Sachs metric describing gravitational radiation from an isolated 
source [|5]||. 

The essential physics underlying the Drummond-Hathrell mechanism for superluminal 
propagation is a violation of the strong equivalence principle. Here, we understand the 
weak equivalence principle (WEP) to be the requirement that spacetime is Riemannian 
and thus has at each point a local inertial frame (LIF). By the strong equivalence principle 
(SEP), we mean the further assumption that the laws of physics are the same in the LIFs 
at each point of spacetime, and take their special relativistic form at the origin of each 
LIF. In particular, this condition states that matter couples to the gravitational field only 
via the connection, with no direct curvature coupling. It is clear that the two forms of 
the equivalence principle have a quite different status. While the WEP is fundamental 
to the structure of general relativity, the SEP appears to be merely an extra dynamical 
assumption (minimal coupling) which may not be essential for the self-consistency of the 
theory. Of course, this is precisely what we are testing by studying superluminal propa- 
gation and causality in a situation where the SEP is relaxed. If the theory is indeed still 
self-consistent, it then becomes an experimental question whether or not SEP-violating 
interactions exist and with what strength. This paper is concerned with the characteris- 
tics of photon propagation in a particularly relevant gravitational field, that created by a 
time-dependent, isolated, radiating source. 

The particular SEP-violating interactions we consider here are given by the effective 
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action 
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As shown by Drummond and Hathrell, this action arises through vacuum polarisation ef- 
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fects at one-loop level in QED in curved spacetime. In this case, a = — ttj — , & = w ~ an< ^ 



c = — ^tt — and the scale m is the electron mass. Further corrections to the effective action 

360 7r 

involving higher derivatives of the field strengths and curvatures also arise in QED and are 
relevant to the question of dispersion and high-frequency propagation, but these will not 
be considered here. Similar effective actions may also be expected to arise generically as 
low-energy approximations to more speculative fundamental theories of quantum gravity. 
However, at this point we do not need to restrict ourselves to any particular mechanism, 
but can instead consider the effective action ( |1 . 1|) as a phenomenological model involving 
a new fundamental scale m to be determined by experiment. This is the approach we will 
adopt in this paper. 

As reviewed in the next section, this action induces curvature-dependent modifications 
to the effective light cone governing photon propagation and leads to the possibility of 
superluminal velocities. Two important generic features were formalised in ref.|Q]. First, 
for Ricci flat spacetimes there exists a polarisation sum rule whereby if one polarisation 
state propagates with velocity less than c, the other must necessarily have a velocity 
greater than c. Second, for black hole spacetimes, even if the photon velocity along some 
trajectory is different from c, it reverts to the standard light cone velocity precisely on 
the event horizon, ensuring that the physical and geometrical horizons coincide. These 
features are most apparent in the following formula for the effective light cone, following 
from eq.(|l.l|): 

k 2 = ~^(b + 2c)T^k»k» + ^-C^ up k»k»a x aP (1.2) 

where k^ is the wave vector (photon momentum) and a x is the polarisation. In this form, 
we have used the Einstein equations (with G = 1) to substitute the energy- momentum ten- 
sor T pu for the Ricci tensor. The first term is polarisation-independent and is proportional 
to the projection of the energy-momentum tensor appearing in the weak energy condition, 
viz. T^k^k 11 > for any null vector k^. A similar modification to k 2 also arises for photon 
propagation in other modified environments such as background electromagnetic fields or 
finite temperature. The second contribution, which is special to gravitational backgrounds, 
produces a polarisation dependent shift in the effective light cone proportional to a partic- 
ular projection of the Weyl tensor (i.e. a particular Newman-Penrose scalar, as explained 
in section 3). The shift is equal and opposite for the two polarisations, one of which is 
therefore necessarily superluminal. 
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A polarisation-dependent shift in the light cone (gravitational birefringence) allowing 
superluminal propagation is therefore the essential signature of SEP-violating interactions 
in curved spacetime electrodynamics. 

If this effect is indeed observable, there are a number of potentially important cases 
which should be considered. The first, already discussed in ref.[Q, is the polarisation 
dependence of the bending of light by a star or other matter distribution whose exterior 
field is described by the Schwarzschild metric. In this case, the angle of deflection becomes 
dependent on the photon polarisation. Superluminal propagation would be especially 
important in the evolution of the early universe, where it could be relevant for the horizon 
problem p|,p|,p|,p!0[] . This probably requires changes in c by many orders of magnitude, 
rather than the perturbative corrections considered here (we are implicitly considering 
( p. . 1|) to be the leading terms in an expansion in 0(R/m 2 )). Nevertheless, ( |1.2|) does predict 
superluminal propagation for FRW spacetimes |]J • Here, the spatial isotropy ensures that 
the second term vanishes (C^\ up = for FRW spacetime), but the first term produces 
a polarisation-independent shift in the light cone which implies superluminal velocities 
provided the weak energy condition p + p > is satisfied. On the other hand, for de Sitter 
spacetime, k 2 remains zero because the spacetime is isotropic. 

Here, we consider a third important case: the propagation of light in the background 
spacetime corresponding to a time-dependent, isolated source emitting gravitational ra- 
diation. We show that there is a polarisation-dependent superluminal effect and express 
this in terms of the appropriate Newman-Penrose scalars characterising the Weyl tensor. 
The asymptotic behaviour depends on the direction of the photons. For ingoing pho- 
tons, the light cone shift is of 0(1 /r) and is given by the amplitude of the gravitational 
waves far from the source; for the more important case of outgoing photons, the shift is 
0(l/r 5 ) and depends on the 'quadrupole aspect' of the Bondi- Sachs metric. These results 
are related to the Peeling Theorem ||, which describes the asymptotic dependence of the 
Newman-Penrose scalars in the Bondi- Sachs metric. 

Before deriving these results, we mention briefly two interesting recent developments 
in the theory of superluminal propagation. First, Albrecht, Barrow and Magueijo |]8|,p|,p!0| 
have studied the implications for the early universe of allowing the fundamental constant c 
to be time-dependent. Although this formalism is very different from that considered here, 
these papers illustrate well the potential importance of superluminal effects in cosmology. 



Second, Drummond |Tl]] has recently proposed a manifestly covariant bi-metric theory of 
gravity in which matter and gravity couple to different vierbeins whose relative orientation 
is determined dynamically through a sigma model action. He then applies this theory 
to the specific problem of 'dark matter', arguing that the phenomena usually attributed 
to the existence of dark matter can alternatively be described by a modification of the 
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fundamental theory of gravity. This is of course a much more radical generalisation of 
conventional gravity than considered in this paper, but this new theory would certainly 
incorporate the type of effects discussed here. 

In the rest of the paper, we review briefly the derivation of the modified light cone 
condition ( |1.2|) in section 2. Then, in section 3, we describe the Bondi-Sachs metric and 
the Peeling Theorem. The results on photon propagation and our general conclusions are 
given in sections 4 and 5. 



2. Superluminal Propagation in Gravitational Fields 

The characteristics of photon propagation following from the effective action ( p. . 1|) are 
most easily described using geometric optics JIJ. In this picture, the electromagnetic field 
is written as the product of a slowly-varying amplitude and a rapidly-varying phase, i.e. 

Fpu = foe" (2.1) 

The wave vector (photon momentum) is defined as = d^d, while the Bianchi identity 
constrains L to have the form 

j 'fj,u = k^ci^ — k u a^ (2-2) 

where the direction of a M specifies the polarisation. For physical polarisations, k^a^ = 0. 

For conventional curved spacetime QED based on the usual Maxwell action, the equa- 
tion of motion is simply 

D^=0 =}► V = (2.3) 

Since this implies 

fcV = (2.4) 

we immediately deduce that k 2 = 0, i.e. k^ 1 is a null vector. Then, from its definition as a 
gradient, we see 

k^D^k" = k^D^kf, = ^D u k 2 = (2.5) 



Light rays (photon trajectories) are defined as the integral curves of the wave vector, i.e. the 

ils 



curves x M (s) where = These curves therefore satisfy 



o = w = ^ + r ^^. (2 . 6) 
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which is the geodesic equation. So in the conventional theory, light rays are null geodesies. 

The effective action ( |1 . 1|) gives rise to a modified equation of motion which, under the 
approximations listed below, implies 



2bR» x kJ x » + 4cR^ Xp kJ x ? 



(2.7) 



Here, we have made the standard geometric optics approximation of neglecting derivatives 
of f^v relative to the derivatives of the phase factor, which produce powers of the mo- 
mentum k^; we have neglected derivatives of the curvature tensor, which would produce 
corrections of 0(X/L), where A is the photon wavelength and L is a typical curvature 
scale; and we have omitted the corrections to the term D^F^ of 0(aR/m 2 ) coming from 
the Ricci scalar in the effective action - these are of 0(A^/L 2 ), where A c is the Compton 
wavelength corresponding to a particle of mass m (the electron in the conventional QED 
derivation of (|1 . 1|) ) which should be neglected as higher order if we view (|1.1[) as the 
leading term in an expansion in powers of curvature. 

Eq.( |2.7| ) can now be rewritten as an equation for the polarisation vector a M , and 
re-expressing in terms of the Weyl tensor we find 

k 2 a u + ( 2b + 4c ) R J^k x a v -k»k u a x ) + ^ C" Xp Wk x a p = (2.8) 
wr m z 

The solutions of this equation describe the characteristics of propagation for a photon of 
momentum and polarisation a?. Contracting with (and assuming spacelike normal- 
isation a M a M = —1), we find the effective light cone 



- k - —Ls^xpK 1 k a a p = (2.9) 



k + 7) Raxk^k KCauXpk^k a a 

m m 

from which (|1.2|) follows immediately. 

It should be noted that all these equations are manifestly local Lorentz invariant. 
On the other hand, the presence of the explicit curvature coupling in the effective action 
means that the equations of motion do not reduce to their special relativistic form at the 
origin of each LIF, and thus that the dynamics is different in the LIFs at different points in 
spacetime. In this sense, these equations violate the strong principle of equivalence. Some 
implications of this for causality have been discussed in refs.p|,|i|,|T2"||. 

At this point, it is illuminating to re- write the effective light cone condition using 
the Newman-Penrose tetrad formalism, and in particular to show the dependence on the 
NP scalars characterising the Weyl tensor. The first step is to choose a null tetrad as 
follows. Let l p be a null vector. Let a M and 6 M be spacelike, transverse vectors and define 
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the complex null vectors m M and m M by m M = 7^(0^ + ib^) and m M = — 

Finally, choose a further null vector n M orthogonal to and m M . These vectors satisfy 
the orthogonality conditions: 

£.m = l.m = n.m = n.rh = (2-10) 

they are null vectors: 

£.£ = n.n = m.m = rh.rh = (2-H) 
and we impose the normalisation conditions: 

t.n = 1 m.m = -1 (2.12) 



The null tetrad is defined by the vierbeins e^, where we define = £ M , = n^, = m A 

A* 

4 



and e 1 ^ = m M . The corresponding metric is 
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(2.13) 



The five complex NP scalar s characterising the Weyl tensor are (following closely the 
notation of ref. [O] ) 



*0 = 


-C abcd £ a m b £ c m d 


= — C1313 


tfl = 


-C abcd tn b £ c m d 


= — C1213 


#2 = 


-C abcd £ a m b m c n d 


= — C1342 


*3 = 


-C abcd £ a n b m c n d 


= — C1242 


V& 4 = 


-C abcd n a m b n c m d 


= — C*2424 



(2.14) 

The symmetries of the Weyl tensor imply several interesting relations amongst its compo- 
nents. Most important for our discussion are the trace-free conditions 

V ad C ab cd = (2.15) 

and cyclicity, e.g. 

C1234 + C1342 + C1423 = (2-16) 
Together, these imply the important identity 

C^ Xp £»m v £ x m p = C1314 = (2.17) 

6 



Components of the Ricci tensor have a similar classification. We define 



$00 = ~R llv in v = ~Ru $22 = -\R^n»n» = -\R 22 

$02 = -^R^m^m" = --#33 $20 = -\R^rh^fh u = -\R^ 

$oi = ~R^m v = -Ir 13 $ 10 = -\R^my = -\R 1A 

$12 = —R^m" = --R23 $21 = -\R^m v = -\R 2A 

$11 = -J(VV + ^ m " m '') = -^(#12 + ^34) 

A = Ya r = ^( R ^ n "- R ^ m ^) = ^(^12-^34) 

(2.18) 

We can now re-express the light cone condition in NP form, using the identity ( |2.17| ) . 
For example, if we choose the (unperturbed) photon momentum in the direction of the null 
vector £ M , i.e. k^ 1 = uj£^, and the transverse polarisation vectors as a M , then the term 
in ( p.9| ) proportional to the Weyl tensor becomes, for the two polarisations respectively, 



±i($ + $*) (2.19) 



The lightcone condition is therefore simply 



depending on the polarisation. It is interesting that this depends on only a single NP 
scalar for each of the Ricci and Weyl tensor contributions. The polarisation sum rule and 
horizon theorem described in ref.[|J are immediate consequences of this form of the light- 
cone condition. Other choices of the photon momentum and polarisation give analogous 
expressions for the modified light cone. 

In the rest of this paper, we apply this result to the special example of the Bondi-Sachs 
metric and show the precise relation to the Peeling Theorem. 
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3. Bondi-Sachs Metric and the Peeling Theorem 



The spacetime describing an isolated, radiating source is given by the Bondi-Sachs 
metric: 

ds 2 = -Wdu 2 - 2e 2p dudr + r 2 h tJ (dx l - U l du)(dx J - U j du) (3.1) 

where 

h ij dx i dx j = \{e 21 + e 2S )d9 2 + 2 sinh( 7 - 5) sinddOdcj) + \{e~ 21 + e~ 2S ) sin 2 6d(j) 2 (3.2) 

The metric is valid in the vicinity of future null infinity I + . The family of hypersurfaces 
u = const are null, i.e. d^ud^u = 0. Their normal vector satisfies 

if, = d^u I 2 = 0, l^Dpt = (3.3) 

The curves with tangent vector £ M are therefore null geodesies; the coordinate r is a radial 
parameter along these rays and is identified as the luminosity distance. 

The six independent functions characterising the metric have the following asymptotic 
expansions near X + for large r: 

£ = -7(4 + 4)4 + 0(1) 



4 X ^ A/ r 

5(7 + a) = ^ + £ + o£) 
5(7-0 = ^ + ^ + 0^) 

U e + i sin eU* = \- (do ^ds) ( sin 2 9(c+ + " ' ' 

sin 9 V sm 9 J \ 



ic 



X 



+ 2(/ + zsin^ + ...)l + 0(l) (3.4) 



where M., c + ( x ), ?+( x )j d 9 ^ are all functions of (u,9,(f)). The form of these expansions 
follows from a careful analysis of the characteristic initial-value problem for the vacuum 
Einstein equations (see refs.||[| and for a textbook review [|T4|). 



The function M.(u,9,4>) may be called [|14| the 'mass aspect'. Its integral over the 
unit 2-sphere, 

M(u) = -3- f dn 2 M(u,9,<t>) (3.5) 

47T J 

represents the mass of the system at T + and is the familiar Bondi mass. Similarly d 6 ^ and 
q+(x) are the 'dipole' and 'quadrupole aspects' respectively. M., d 9 ^ and q+(x) satisfy 
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dynamical equations derived from the Einstein equations for their u— derivatives, involving 
also the remaining functions c+f x y 

<9 u c_|_(x) are functions which must be given as initial data and are specified on X + . 
They are the Bondi 'news functions'. An especially important result relates d u M.{u) to 
the news functions: 



d u M(u) = ~J dn 2 ({d u c + ) 2 + (d u c x ) 2 ) 



(3.6) 



This states that the Bondi mass is reduced if the news function is non-zero, corresponding 
to the fact that the system loses mass if and only if it is radiating. 

Finally, as we discuss later, the second derivatives d 2 c + ( x ) can be identified as the 
amplitude of the gravitational waves in a weak-field limit. Notice also that 7 ±5, and hence 
c + ( x ) and q+( x ), correspond to the two independent gravitational wave polarisations. 

We return to the interpretation of c + ( x ) in section 4 after introducing the Peeling 
Theorem, discussed by Sachs in ||. This gives the leading asymptotic behaviour in 1/r 
of the set of Newman-Penrose scalars ^0, • • ■, ^4 characterising the Weyl tensor in the 
Bondi-Sachs spacetime. Although we do not need to exploit it here, the • • - ; ^4 are 
intimately related to the Petrov classification, which classifies the Weyl tensor according 
to the degeneracy of its principal null vectors. 

For our purposes, we simply need the following result (note that since we only need 
the precise results for ^0 and ^4 in what follows, we have only written schematic forms 
for the others): 



^4 
^3 

^0 

where we have set £, 



1 



1 



2C> 



sin 2 6 
M+[d e -\ 



sin 9 

i 

sin 9 



sin 2 9d u (c + — ic> 



2 1 2 
c + + c x 



d + i sin 
6(q + - iq 



d^u and chosen the transverse vector 



(3.7) 



^(8^9 + 1811198^ 



The essence of the Peeling Theorem is the correlation between the leading order in 
1/r and the type of the NP scalar. Notice also that the leading coefficients as we pass 
from \l/4 to ^0 involve respectively d 2 c, d u c, M., d and q, with the higher moment aspects 
being associated with successively higher powers of 1/r. 
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4. Photon Propagation in the Bondi-Sachs Spacetime 



We are now ready to combine the general results for modified photon propagation in 
section 2 with the special features of the Bondi-Sachs gravitational radiation spacetime. 
Consider first the case of photons with momentum = ujI^. This corresponds to motion 
radially outwards from the gravitationally radiating source. (To confirm this, note that 
the equipotential surfaces for outgoing waves are $(u) = const, so using the geometrical 
optics analysis above, the corresponding rays have tangent vector /c M = d^fi, so we can 
identify the directions of /c M and Choose the transverse polarisation vectors a M , 6 M to 
lie in the 9 and <f> directions respectively, so that 

a" = —^(m** + m M ) 
V2 

b» = l -{m^-rh^) (4.1) 

V 2 

with m M as in section 3. In this case, according to eq.( [2.20|) the light-cone shift for the 
two polarisations is 

k 2 = ±^(¥ + *5) (4.2) 



which for the Bondi-Sachs metric is 



e = ± ^«± to) 

corresponding to a velocity shift Sv proportional to ±q+/r 5 . As usual for Ricci-flat space- 
times, the two transverse polarisations have equal and opposite velocity shifts, so one is 
always superluminal. For outgoing photons, therefore, we do find a velocity shift, but 
it is very weak, falling off as 1/r 5 , and is governed by the quadrupole aspect q + of the 
gravitational field. 

It is interesting also to look at the photons with polarisation vectors rotated through 
45°. In this case, we choose 

1 i 
a M = -(m M + m M ) (m M - m M ) 

W = --(m^ + rh fl ) - -(m M - m") (4.4) 



A calculation along the lines of eq.( 2.19 ) now gives 

, 1 . 4:CU 2 



m 2 



% 
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For the Bondi-Sachs metric this is 

e = (4.6) 

As expected, the photons with polarisations rotated through 45° are influenced by the x 
polarisation of the gravitational radiation, compared with the original choice aligned with 
the + polarisation. 

A significantly larger effect is obtained if we consider incoming photons, moving ra- 
dially towards the source of gravitational radiation. In this case, we have /c M = un^. For 
the initial choice of photons polarised in the 6 or directions, we find 

k 2 = ^~ (* 4 + n) (4.7) 

Tfl 

while for the polarisations rotated through 45° we have 

e = ±^(* 4 -**) (4.8) 



In Bondi-Sachs, this means 



and 



, 9 8cuj 2 1 _ 9 . , 

k 2 = ±— =- -d 2 u c+ (4.9) 



m 2 r 



k 2 = ±^-d 2 u c x (4.10) 



m 2 r 



respectively. In this case, the superluminal velocity shifts are of 0(l/r), i.e. 5v is propor- 
tional to <9 2 c + ( x )/r depending on whether the photon polarisation is aligned with the + 
or x gravitational radiation polarisation. 

We now see clearly the relation of superluminal photon propagation to the Peeling 
Theorem. Depending on the direction (and polarisation) of the photons, the shift in the 
light cone is proportional to one of the NP scalars characterising the Weyl tensor. The 
Peeling Theorem specifies the leading order in 1/r of each of these types in the vicinity of 
X + , which translates immediately into a result giving the 1 /r-dependence of the photon 
velocity shifts. 

These results of course hold for the full gravitational radiation field described by the 
Bondi-Sachs metric. It is interesting at this point to compare them with those obtained 
previously for weak-field gravitational radiation in the linearised approximation. To see 
this relation, consider the following metric describing gravitational plane waves in the 
linearised, weak-field limit where the metric perturbation from Minkowski spacetime 
is chosen in transverse, traceless gauge: 

ds 2 = -dudv + (l - h + (u))dx 2 + (l + h + (u))dy 2 - 2h x (u)dxdy (4.11) 
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h + and h x of course correspond to the + and x polarised gravitational waves, and d 2 h + ^ 
represent their amplitudes ||14|| . The relevant components of the Weyl tensor are: 



C U XUX C U y U y o^!l^+ 



1 

2' 
1 



C uxuy = --d 2 h x (4.12) 

Linearising the Bondi-Sachs metric in 7 and 5 and discarding the functions W, f3 and 
U % , we find 

ds 2 = -dudv + r 2 ((1 + 7 + 5)d6 2 + 2(7 - S) sinOddd^ + (1 - 7 - S) sin 2 ## 2 ) (4.13) 

which has therefore reduced to the weak-field gravitational wave metric with 

h+ = -( 1 + S) = -2^± h x =-(7 -5) = -2^ (4.14) 

confirming the identification already assumed above that c + ( x ) correspond to the two 
independent gravitational wave polarisations. 



Returning to the metric ( f4.11| ), the light-cone shift for photons travelling in the op- 
posite direction to the gravitational waves is 

k 2 = ±^C UXUX = ±^d 2 u h + (4.15) 
for x, y polarised photons, and 

k 2 = ±^C uxuy = ±^-dlh x (4.16) 

for 45° rotated photons. With the weak- field identifications ( |4.14|) , we recover ( |4.9|) , ( [4.10|) . 

On the other hand, for photons travelling in the same direction as the gravitational 
waves, the light-cone shifts are proportional to the C vxvx , C vyvy and C vxvy components of 
the Weyl tensor, which vanish. For both polarisations, therefore, k 2 = 0. 

To summarise, for weak-field gravitational plane waves, the effect on photon velocity is 
as follows. Photons travelling in the same direction as the gravitational waves feel no effect, 
whereas photons travelling in the opposite direction experience velocity shifts proportional 
to the gravitational wave amplitudes <9 2 /i_|_( x ), depending on the alignment of the photon 
and gravitational wave polarisations. 

In the full Bondi-Sachs radiation metric, a similar situation holds, except that here 
the photons travelling in the same direction as the gravitational radiation also experience 
a velocity shift proportional to q+( x )i though only of 0(l/r 5 ), while those travelling in the 
opposite direction experience a much larger 0(l/r) effect proportional to <9 2 c_|_( x ). 
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5. Conclusion 

If the strong equivalence principle is violated, photons do not necessarily propagate 
along the null geodesies of the background curved spacetime. The physical light cone 
is shifted with respect to the geometrical one. Here, we have discussed this effect in 
terms of an effective field theory containing explicit SEP-violating interactions. Such 
terms have been shown to arise even in conventional QED in curved spacetime through 
vacuum polarisation effects and are expected to appear generically in the low-energy limit 
of theories of quantum gravity. 

Developing earlier work on black hole spacetimes, we have considered the special case 
of the gravitational radiation metric introduced by Bondi, van der Burg and Metzner || 
and by Sachs . This has allowed us to generalise previous results on photon propagation 
in weak-field gravitational wave backgrounds obtained by Drummond and Hathrell ffl. 

Our principal results are that the velocity shifts are either superluminal or subluminal 
depending on which of the transverse photon polarisations is considered. For photons trav- 
elling inward towards the source of gravitational radiation, the velocity shifts are asymp- 
totically of 0(l/r) and depend on the functions <9^c_|_( x ), where d u c + ^ are the Bondi 
news functions. The light-cone shift depends on the relative alignment of the photon and 
gravitational radiation polarisations. For photons moving outwards along the direction of 
the gravitational radiation, we still find a non-vanishing velocity shift (in contrast to the 
case of weak- field gravitational waves) controlled by the quadrupole aspect <?+( x ) of the 
gravitational radiation, but it is only of 0(l/r 5 ). 

The implications of these results for the question of causality in the presence of su- 
perluminal propagation are left for future work. 
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